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Abstract 

The thermal model is commonly used in two different ways for the descrip- 
tion of hadron production in ultra-relativistic heavy ion collision. One is the 
application of the thermal model to Air integrated data and the other is the 
thermal description of central dN/dy ratios. While the first method implicitly 
assumes global equilibrium the other scenario assumes Bjorken scaling within 
the investigated rapidity range. Both assumptions are only approximations 
for real physical collision systems. We study the impact of both approxi- 
mations for the extraction of thermal parameters on the exemplary case of 
S+S collisions at SPS energies. The particle distributions are modeled by a 
hydrodynamical description of the relevant collision system. 
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I. INTRODUCTION 



One goal of the studies of high energy collisions is the understanding of the inclusive 
hadron production [p]]. A relevant question is whether the hadron production is dominated 
by phase space and statistical laws or by dynamical constraints of the underlying theory of 
the strong interaction. There is strong evidence that even in elementary collisions like e + e~ 
or pp at sufficiently high energies the dominant factor is a statistical filling of phase space 

01- 

In heavy ion collisions one may be even closer to thermodynamic behaviour due to 
secondary interactions. Therefore a big effort is going on to study thermal behaviour in 
these collisions by microscopic models as well as to classify directly the experimental 
final hadronic state by a thermal model PHTB"||. While within the thermal model applied 



at low incident energies the whole system may be regarded as one fireball in approximate 
global thermal and chemical equilibrium this is not anymore justified in ultra-relativistic 
heavy ion collisions. The incoming nucleons are only partly decelerated and the leading 
protons cannot equilibrate with the mesons in the center. However, the assumption of local 
thermal equilibrium might still be a valid assumption in these reactions. 

If there is only local thermal equilibrium a fundamental problem arises in analyzing 
particle yields directly from experiment. A measured hadron does not tell us from which 
spatial region it comes from. In principle a model is needed which provides the spatial 
information about the particle production, like the hydrodynamic model or so called event 
generators. On the other side one may try to analyze chemical equilibration of particle 
production in ultra-relativistic heavy ion collisions directly from the measured particles 
without a detailed model about spatial differences in thermodynamic quantities. This is 
done in two different ways where both methods make necessarily some compromise. 

The first method is to restrict oneself to a limited region in momentum space (e.g. 
around midrapidity) and to assume that the considered region in rapidity corresponds to a 
homogeneous spatial region in thermal equilibrium, i.e. all thermal parameters are constant 
over that region [fTr-T^j. This procedure has the disadvantage that there is no one-to-one 



correspondence of spatial regions to regions in momentum (rapidity) space. One local cell 
contributes in general to different rapidity regions and even differently for particles of differ- 
ent mass. However, in the limit of infinite collision energy we have Bjorken scaling jnj which 
assures that along longitudinal proper time contours no spatial gradients of local variables 
are present. (There may still be gradients in transverse direction which we assume to be 
small). We recover global equilibrium again. Note that we use here and in the following the 
notion "global equilibrium" if spatial regions have the same values of all intensive thermo- 
dynamic variables even if there is no causal connection between these regions. In the scaling 
limit a thermal and chemical analysis of particle spectra directly from experiment is possible 
again and in the case of RHIC and LHC this may be the only reasonable way to do it. 

At SPS energies and below the rapidity spectra don't show Bjorken scaling. Therefore a 
chemical analysis in a limited rapidity region might be questionable P0 |. We will investigate 



here how well the Bjorken scaling assumption in connection with a chemical analysis works 
using a hydrodynamic model at SPS energies. 

The mentioned problem of rapidity cuts leads to the second method for a chemical 



analysis, i.e. a global thermal model to An integrated data [16-18|. This method implicitly 
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assumes global chemical equilibrium in the above mentioned sense. This assumption is less 
and less justified the smaller the colliding nuclei are and the higher the collision energy is. 
In S+S collisions at CERN-SPS a clear deviation from global equilibrium between pions and 



net protons is seen in that rapidity spectra pl |. In our study here we first like to address 



the question about the error of a global thermal fit to Air integrated data in the case of only 
local equilibrium. 



II. LOCAL AND GLOBAL EQUILIBRIUM 

Let us assume that particles decouple (freeze-out) from spatial regions which build up 
a continuous 3-dim freeze-out surface embedded in the 4-dim space-time. Along this 
surface we assume further to have local thermal and chemical equilibrium, but the thermal 
parameters may vary along this surface. The local particle densities Pi(x) are then given 
by the local thermal parameters, temperature T(x), baryon chemical potential Pb(%), and 
strangeness chemical potential (J>s(x): 

T(r) 00 m 2 im- 

Pl {x) = (2Jj + 1) -AJ- J2( T iy +1 expbGM*)^ + l*(x)Si)] -?- MlA) , (i) 

Z7r j=l J 1 \ x ) 

where J, is the spin degeneracy, -B, the baryon number, Si the strangeness, m ; the mass 
of particle % and the — (+) sign is for fermions (bosons), respectively. We neglect possible 
suppression factors like strangeness suppression %(x) etc which are needed if only relative 



chemical equilibrium is present |]22| . We also assume that strangeness is conserved locally, 
i.e. the net strangeness density is zero everywhere. Thus fis( x ) can be expressed locally as 
a function of T{x) and Hb(x). The multiplicity of a particle species i is calculated by 



iVf cal = J d^(x)if[T(x),/iB(x),/is(x)], (2) 

E(x) 



where jf is the particle current of species i. We use the definition of Eckart j£3| for the four 
velocity and we can therefore decompose the current (x) = pi(x)-u M (x) into the product 
of the local rest frame density Pi(x) and the four velocity u^{x). 

In the case of a global thermal fit one assumes that all thermal parameters are constant 
and thus the multiplicities are given by 

Nf lobal (T,// B ,A*s) = / d 3 cr At (x) pi(T,/x B , ^s) u^{x) = Pi(T,// B , A*s) V , (3) 

S(x) 

where V — / s d 3 <r A1 (a;) u^{x) is the Lorentz invariant comoving eigen volume. We define the 
freeze-out average of a thermal parameter O by 

{Q) = ; s d 3 ^(x) vHx) qjx) _ (4) 

Next we define the global fit thermal parameters &t as the parameters which give the 
minimum x to experimental data or in our study case to the iVj ocal by 
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X 2 (^ fit , T fi \ ^) = £ ^ 1 ' = Min X 2 (^, T, „ B ) . (5) 



For small variations of the freeze-out parameters along the freeze-out surface we have 
the following approximate relation: 

A^(dV^)) » iVf lobal (7*,^,/^) « iVf lobal ((T>, (^b), (/*)) . (6) 

However, in general all three quantities of Eq. (||) are different. Here we like to study the 
validity of Eq. (|6|) for a realistic example in order to get a feeling about the goodness of a 
global fit to particle yields which arise from a system which is only in local but not in global 
equilibrium. 



III. THE HYDRODYNAMICAL STUDY CASE OF S+S COLLISIONS 

The minimal deviation from global equilibrium forced by the experimental spectra is that 
the baryon density differs locally in space. So far there is no convincing evidence of a local 
change in temperature and therefore it is usually assumed to be constant (but see also p5fl ). 
A realistic freeze-out surface with constant T but varying baryon density and thus varying 
lib and [j,s is naturally provided by a hydrodynamical simulation p6| . Hydrodynamics by 
definition assumes local thermal and chemical equilibrium. 

We take as an example for our studies S+S collisions at CERN-SPS since this collision 
system shows the largest discrepancies in the proton and pion rapidity distribution and 
therefore deviations from global equilibrium. We perform a hydrodynamical simulation of 
S+S collisions with the same initial conditions and in the same way as described in reference 
f2?|| . The only difference is that the freeze-out hypersurface is defined now on the contour 
of constant temperature T = 140 MeV. The temperature and the chemical potentials follow 
from the local energy density and baryon density with the help of the used equation of state 
which was defined in |28| and labeled as EOS A. This equation of state contains very few 
hadronic resonances. In order to have a typical resonance spectrum for a chemical analysis we 
use the resonance spectrum up to a mass of 1.7 GeV for calculating the particle spectra. This 
introduces a small inconsistency since the equation of state in the hydrodynamical evolution 
is different from the equation of state used for particle spectra. However, at the low freeze- 
out temperature of 140 MeV the higher resonance states are of minor importance. We show 
that the calculated freeze-out hypersurface is still compatible with the higher number of 
resonance states by comparing the resulting spectra with experiments. In Figure [l] we show 
as solid lines the h~ and net proton rapidity and transverse mass distributions. The spectra 
are calculated as described in p8| using the description of Cooper and Frye f29| . We see 



that the calculated spectra are still in reasonable agreement with the data despite the use 
of the larger resonance input. 

We now have a model system which is clearly out of global chemical equilibrium. In 
order to show the deviations from global equilibrium, we plot in Figure |2| the distribution of 
sub- volumes dV/d/ZB as function of /z B as they result from our hydrodynamical simulation. 
The width in [i B is of order 100 MeV around the average of (lib) = 193 MeV. We see a large 
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spread in \ib indicating that there is a large deviation from global chemical equilibrium with 
respect to baryon number. 

The resulting spectra are integrated over 4tt and the multiplicities Nl ocal are shown in 
Table |I[ In order to perform a \ 2 fit we give these yields a relative error as typical for 
experiment [5(J shown in brackets in Table |. 

In a first attempt to describe the local yields one may take a global equilibrated thermal 
model with parameters resulting from averaging over the freeze-out hypersurface. We first 
take as chemical potential for the global description the averaged values (/xb) and (/is)- The 
result is given in the third column of Table Q. All yields are underpredicted. Mesons come 
out right but baryons and anti-baryons yields are too small. The reason is that baryon yields 
are proportional to the baryon fugacity Ab = exp(/iB/T) and (/ib) < TIu(Ab)- In the case 
of the anti-baryons we have exp(— ((ib)/T) < (Ag 1 ). 

Using the average chemical potential leads in a global model to a reduction in the total 
baryon number. In order to avoid this problem one may use (Ab) and (As) as parameters 
in the global model. Since the baryon yields are proportional to Ab in Boltzmann approx- 
imation the local and global numbers for non-strange baryons are the same up to minor 
corrections due to resonance decays and Fermi statistic. The result of such a calculation is 
shown in the fourth column of Table |. Such a scenario, however, leads to large discrepancies 
for the anti-baryons. Therefore such a description is not satisfactory, either. 

Now we perform a fit to the local yields with a global thermal model. We take as fit 
parameters the volume V, temperature T, and the baryon chemical potential /jb- A*s is 
determined by the requirement of strangeness neutrality and not used as a fit parameter. 
The result of the fit is shown in Table |, too. We recover in this fit nearly the input 
temperature and get a ji^ = 205.7 MeV which is between the average (/ib) = 193.2 and the 
/ib resulting from the average (Ab), £*b — 218.4 MeV. The deviations of individual yields of 
the global fit from the local integrated ones are small. The average deviation is of order 4%. 
The largest deviations are of order 10% for the anti-nucleons and the fl. 

We conclude that the performance of a global thermal fit to Att integrated data is fine 
because the deviations in temperature and volume from the exact numbers are small in the 
studied case of S+S collisions at SPS energies. We expect that going to larger nuclei and 
to smaller energies the amount of stopping increases and therefore the assumption of global 
equilibrium for extracting thermal parameters is even more reliable. 



IV. RAPIDITY CUTS 

Next we study the influence of cuts in rapidity on the extraction of thermal parameters. 
For all particles in Table | we integrate the corresponding spectra of the hydrodynamical 
simulation only over a finite interval in rapidity similar to our studies in |20| . The resulting 
particle yields are fitted in the same way as done before in case of Air yields. This means 
that we assume Bjorken scaling in the sense that the multiplicities of particles in a finite 
rapidity range are still given by Eqs. (0|1|). The resulting thermal parameters are shown in 
Table |I|. Before discussing the result of that exercise we construct two hypothetical cases for 
particle production in order to compare with. First we take the result of the hydrodynamical 
simulation and give every fluid cell on the freeze-out surface by hand a constant /iB = 
140 MeV x In (Ab) = 218.4 MeV. Then we have a system in global equilibrium, but still 
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exhibiting the same flow in longitudinal and transverse direction as in the hydrodynamical 
simulation. Some of the resulting spectra are shown in Figure |l|. The spectra are integrated 
over finite rapidity intervals and fitted with the global thermal model as done before. The 
resulting thermal parameters are shown in Table [II. Table |IV| shows the result of the 
same procedure applied to the spectra of a static fireball with the same volume as the 
hydrodynamical simulation. The spectra are also shown in Figure [l]. 

In case of the static fireball with no flow, we see the largest influence of the cuts. The 
extracted temperature and /xb change considerably. Also the \ 2 increases drastically going 
to smaller rapidity windows. We conclude that a thermal analysis in a limited rapidity 
interval for a static fireball is unreasonable if the analyzed rapidity interval is smaller than 
the thermal width of the lightest particle. 

In the case of global equilibrium with longitudinal flow the picture changes. Even though 
we don't see Bjorken scaling in the rapidity spectra of Figure |I| the extracted thermal 
parameters are rather constant and the quality of the fit stays acceptable even for the 
smallest rapidity window (see Table pip . We see a tendency of increasing temperature and 
/iB with decreasing rapidity interval. This increase is artificially induced by the rapidity 
cuts but much weaker than in the case of the static fireball. Since at AGS energies and 
especially at SIS energies we expect less longitudinal flow, the artificial increase of the fitted 
temperature due to rapidity cuts around midrapidity may be larger. 

In the hydrodynamic case of Table |I| we see a drastic decrease of /is due to the baryon 
hole at midrapidity. The temperature, however, shows a very similar, only very small in- 
crease as in the case of global equilibrium with flow and may therefore be attributed to an 
artificial increase due to rapidity cuts. The quality of the fit is rather independent of the 
cut. In the same exercise [20] with the yields from RQMD |3(J we saw a larger increase of 
X 2 with decreasing rapidity window and a larger change of thermal parameters e.g. tem- 
perature. This is due to the fact that RQMD yields are not in perfect local thermal and 
chemical equilibrium as it is assumed here. Especially the strange hadron production is 
quite different in the central region compared to the fragmentation regions. In other words, 
in the study of rapidity cuts within the RQMD model purely kinematic bias on thermal 
parameters cannot be separated from the impact of different physics in central regions com- 
pared to fragmentation regions. This is different from our study here, where any changes in 
the thermal parameters of tables [TT1] and [TV] are artificial changes due to improper kinematic 
cuts. 

We summarize that a thermal fit to yields or ratios in a limited rapidity region gives 
reasonable results as long as there is large enough longitudinal flow. 



V. CONCLUSION 

We have studied the impact of two approximations often used in the thermal analysis of 
experimental particle ratios and yields. We have done this studies for the exemplary case of 
S+S collisions at 200 A GeV. First we have shown that in the case of only local equilibrium a 
thermal description of the 47r integrated data by a global thermal model leads to deviations. 
The reason is that in general the average particle density (pi) is different from the density 
resulting from average parameters, i.e. 
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{ Pi [T(x) , /x B (x) , /xs (x)] ) ^ Pl ( (T) , (/! B > , (/i S ) ) 



(7) 



The deviations depend on how the average is taken and can be as large as 20%. 

Minimal deviations are achieved by performing a x 2 -fit to the 4n integrated yields from a 
hydrodynamic simulation with local chemical equilibrium. Such fits reproduce the constant 
input temperature up to a few percent and lead to a /xb which is of the order of the average 
(/xb)- The deviations are generally small (up to 10%). The quality of a global thermal 
model in case of only local equilibrium in other cases than studied here, e.g. local variations 
of temperature, larger variations in the local baryon density, etc. have to be investigated 
individually. However, we think that the result of a reasonable description of yields from a 
only local equilibrated system by a global model will to a large extend remain valid. 

We also studied the influence of rapidity cuts on the extraction of thermal model param- 
eters. We explicitly showed that in a system without longitudinal flow the rapidity cuts lead 
to serious problems. In the case of S+S collisions at SPS, however, the cuts don't spoil the 
extraction of the temperature but lead to smaller central /xb as it is expected from the dip 
in proton rapidity spectra. We conclude that at SPS energies already enough longitudinal 
flow is present to justify the Bjorken scaling assumption, in which case the fitted thermal 
parameters are independent of rapidity cuts. 

The decision which of both methods should be used for a chemical analysis depends on 
the amount of longitudinal flow in the system. For low energies and large systems a 4tt 
analysis is recommended while for small systems and high energies the analysis of dN/dy 
around midrapidity should be done. We have shown that for S+S collision at SPS energies 
both methods give reasonable results. For RHIC and LHC an analysis in the central dN/dy 
is recommended while for lower energies, especially at GSI we strongly recommend to analyze 
the Air integrated data. 
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TABLES 



TABLE I. 4-7T particle yields: hydrodynamical simulation of S+S collisions with local varying 
baryon density (second column); global thermal model using the average (//b) an d (ms) (third 
column) or (Ab) and (As) (fourth column) of the hydrodynamical simulation; global thermal fit 
(fifth column) to the yields of the second column. The errors in brackets of column 2 are used for 
the x 2 n t and are motivated by typical experimental errors j^J. The numbers in square brackets 
give the percentage of deviation from mean numbers of the local yields of column 2. 
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Zj 




fl 1QQ [ft QO^l 


117 [91 


u /oj 


0.141 




[-5.1%] 


H° 


0.659 (±0.066) 


0.630 [-4.4%] 


0.742 [12 


5%] 


0.710 




[7.7%] 


H° 


0.097 (±0.014) 


0.094 [-3.5%] 


0.080 [-18 


1%] 


0.097 




[-0.1%] 




0.644 (±0.064) 


0.615 [-4.4%] 


0.725 [12 


5%] 


0.694 




[7.8%] 




0.095 (±0.014) 


0.092 [-3.5%] 


0.078 [-18 


1%) 


0.095 




[-0.0%] 


n- 


0.073 (±0.017) 


0.072 [-1.6%] 


0.084 [14 


9%] 


0.082 




[12.3%] 


Q- 


0.017 (±0.006) 


0.016 [-1.5%] 


0.014 [-15 


7%] 


0.018 




[6.0%] 


T (MeV) 


(140.0) 


140.0 


140.0 




141.8 




(±1.2) 


V (fm 3 ) 


1220 


1220 


1220 




1134 




(±71) 


MB (MeV) 


(193.2) 


193.2 


218.4 




205.7 




(±5.3) 


Ms (MeV) 


(29.9) 


29.9 


31.0 






31.8 




A B 


(4.758) 


3.975 


4.758 






4.26 




As 


(1.248) 


1.238 


1.248 






1.255 


x 2 




19.89/22 


15.35/22 




3.97/22 
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TABLE II. Thermal fit to particle yields in various rapidity intervals resulting from a hy- 
drodynamical simulation of S+S collisions at constant freeze-out temperature of T = 140 MeV 
corresponding to the solid line in Figure |l|. Only the resulting thermal parameters are shown. 

parameter y <0.5 y <1.0 y <1.5 y <2.0 4tt 

T(MeV) 143.4±1.2 143.5±1.2 144.0±1.2 144.0±1.2 141.8±1.2 

V(fm 3 ) 316±19 594±36 793±49 932±59 1134±71 

/i B (MeV) 154.8±5.0 165.0±5.1 186.2±5.2 205.3±5.3 205.7±5.3 

X 2 3.40/22 2.45/22 2.64/22 3.85/22 3.97/22 



TABLE III. Thermal fit to particle yields in various rapidity intervals resulting from a hydro- 
dynamical simulation of S+S collisions at constant freeze-out temperature of T = 140 MeV and 
taking artificially a constant /xb = 218.4 MeV at freeze-out corresponding to the dotted line in 
Figure |l|. Only the resulting thermal parameters are shown. 



parameter 



y <0.5 



y <1.0 



y <1.5 



y<2.0 



4tt 



T(MeV) 
V(fm 3 ) 
/i B (MeV) 



X 



143.4±1.2 
314±20 

226.9±5.2 
4.37/22 



143.4±1.2 
593±37 

226.8±5.2 
3.53/22 



143.2±1.2 
819±51 

226.2±5.2 
3.07/22 



142.2±1.2 
1006±61 

223.9±5.1 
2.09/22 



140.0±1.1 
1224±73 
218.5±5.0 
0.050/22 1 



1 Here \ 2 should be exactly zero. The finite value is due to the numerical uncertainty resulting from the 
integration of discretized momentum spectra. 



TABLE IV. Thermal fit to particle yields in various rapidity intervals resulting from a static 
fireball with T = 140 MeV, = 218.4 MeV and V = 1220 fm 3 corresponding to the dashed line 
in Figure |l[ Only the resulting thermal parameters are shown. 



parameter 
T(MeV) 
V(fm 3 ) 

MB (MeV) 

X 2 



y <0.5 
155.1±1.5 

343±26 
259.0±6.1 
61.86/22 



V <l-0 
144.7±1.2 

860±54 
230.7±5.3 
13.44/22 



y <i-5 

140.8±1.1 
1154±69 

220.6±5.0 
1.16/22 



V <2-0 
140.0±1.1 

1225±73 
218.7±5.0 

0.092/22 



4tt 
139.8±1.1 
1240±74 
218.4±5.0 
0.046/22 2 



2 Here x 2 should be exactly zero. The finite value is due to the numerical uncertainty resulting from the 
integration of discretized momentum spectra. 
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FIGURES 



FIG. 1. Comparison of particle spectra of a hydrodynamical simulation with freeze-out at 
constant temperature T = 140 MeV (solid line) with experimental data of NA35 [24]. The dotted 
line results if all fluid cells at freeze-out is given by hand a constant [ab( x ) = 218.4 MeV and 
corresponds therefore to a scenario of global equilibrium. The dashed line correspond to the 
spectra of a static fireball with the same volume as the hydrodynamical simulations as well as 
T = 140 MeV and /i B = 218.4 MeV. 
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FIG. 2. Distribution of sub-volumes V having a given freeze-out baryon chemical potential 
81 1 1 1 1 1 
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